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ABSTRACT 

Using the general form of the static energy solutions to the Dirac equation 
with a magnetic field, we calculate a general self-energy matrix in the 
Furry-picture. In the limit of high temperatures, but even higher magnetic 
fields, a self-consistent dispersion relation is solved. In contrast to the 
high temperature limit, this merely results in a small mass shift. The 
electron anomalous magnetic moment is calculated. The contribution 
from thermal fermions is found to be different from the corresponding 
contribution using perturbation theory and plane-wave external states. In 
the low temperature limit the self-energy is shown to exhibit de Haas-van 
Alphen oscillations. In the limit of low temperatures and high densities, 
the self-energy becomes very large. 
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1 Introduction 



Several compact stellar objects (e.g. neutron stars, white dwarfs and red giants) and 
cosmological models are characterized by high temperatures and/or densities, and 
also the presence of a strong background magnetic field. It is thus of great interest 
to compute the properties of electrons in such extreme environments. Earlier the 
effective potential has been considered^. Here we shall deal with the electron self- 
energy, important for kinematical issues. Apart from the general self-energy matrix 
in the Furry-picture, and the treatment of the de Haas-van Alphen oscillations in the 
self-energy, this work is merely a shorter version of results published elsewhere^, to 
which we refer for further details and references. 



2 Furry— Picture Propagator 

The Dirac equation for an electron (charge -e) in an external field, with the vector 
potential A'^, reads 

{i^ + e4 - m)¥^\x) = , (1) 

where k denotes a complete set of quantum numbers necessary to specify the solutions, 
and the superscripts -|- and — are referring to positive and negative energy solutions, 
respectively. Using a complete set of static energy solutions we may represent the 
second quantized electron field in the Furry-picture^ 

vI/(x,t)=^[6.M/W(x,t)+4vl/(-)(x,t)] , (2) 

K 

where {b, d) are the standard annihilation operators for particles and anti-particles. 
The fermion propagator, including the effects of some distribution of particles, can 
then be constructed explicitly as the expectation value0 

z5(x',x) = (T [*(x')"*"(x)l) . (3) 



In the case of a static uniform magnetic field in the negative z-direction we may 
choose Ay = {0,0, Bx,0). The solutions to the Dirac equation must then be of the 
formi 

= ^ exp[±i{-Ent + PyV + PzZ)]Vn,py {X)u^,n,py,p, , (4) 

where ^ = ±1 is a polarization index, n = 0, 1,2, . . ., and Py,Pz G R. The energy 
spectrum is given by 



En{B, p,) = Jm? +pI + 2eBn . (5) 



In the chiral representation of the 7-matrices we have 

Vn,pyix) = diag[In,pyix),In-l,pyix),In,pyix),In-l,pyix)] , (6) 

where = (eB /ny/^ exp[-eB /2 (x - py/eBY]l/V^. Hn[V2eB{x - Py/eB)]. 

Here Hn is a Hermite polynomial, and we define I-i,py = 0. The spinor U(^^n,py,p:, is 
independent of x'^. There is a twofold degeneracy corresponding to ( = ±1 in all but 
the lowest Landau level n = 0. The electron propagator is with this choice of gauge 

oo „ 

S{x"',x'') = Y. d:ko(tky(tk, Sn,ko,ky,kAx' ,x) 

n=0 

Kq — — — leBn + le 
X exp[-iko{t' - t) + iky{y' - y) + ik,{z' - z)] , (7) 

where (f"/c = d'^k/{2TT)'^, and /f(^o) = 0(^o)/f(^o) + ©(— ^o)/f (^o)- In the case of 
thermal equilibrium the one-particle fermionic distribution is /^(^o) = l/^e'^^^^^^^ + 
1), where (3 is the inverse temperature and fi is the chemical potential determined 
by the charge density of electrons and positrons of the system. Not specifying the 
representation of the 7-matrices we may writei 

Sn,ko,ky,kA^'i^) = ihlO - hlz + m)[a+In,kyix')In,kyix) + Cr-In-l,kyix')In-l,kyix)] 

-iV2eBn['y+In,ky{x')In-i,ky{x) - 'j-In-i,ky{x')In,kyix)] , (8) 
where 7± = (7^ ± ijy)/2, and a± = (1 ± o-z)/2 ((8)lwhen necessary). 



3 The General Self-energy Matrix 



There are two possible contributions to the one-loop electron self-energy. The tadpole 
contribution is proportional to the total electric charge and current in the medium^ 
and thus vanishing in a neutral environment. Therefore we are left only with the IPI 
self-energy 

— i'E{x',x) = {—ieY'y^iD'^'^{x' — x)iS{x',x)'yi, . (9) 
The photon propagator may be written 



iD^'^(x) 



d_ 

dq^ 



—I 



q^ + is 



^rr6{q')fB{qo) 



(10) 



where the photon distribution function in the case of thermal equilibrium is fsiQa) 
l/(exp[/3|go|] — 1)- The effective Dirac equation may then be written 



'x){t^ + e4-m)^^^l^^^p 



[X 



dtdydzU(;,N,py,p,^N{B,Pz]Py)Ui;',N,py,p, , (11) 



1 

n! 

where we have used the general form of the electron wave- functions in Eq.(^ to define 
a general self-energy matrix S^v. The space-time integral (i.e. energy — momentum 
conservation) imphes that both sides of Eq.(|llD are diagonal in n,py and pz- On- 
shell, i.e. with tree-level wave-functions, the effective Dirac equation { ^A\) results in 
an energy shift E^]l{B,Pz) = En{B,pz) + AEi^^n{B,Pz), where 

iPz) — '^C,n,Py,Pz'^"'i'^'P^''Py)'^(','n;Py,Pz ' 

(12) 

that actually is independent of Py. The on-shell self-energy has also been shown@ to 
be independent of the gauge- fixing parameter ^. If not otherwise stated we shall use 
the Feynman gauge ^ = 0. 



4 Strong Field Limit 

In favorable conditions the gauge-fixing dependence of the self-energy may be ne- 
glected also off-shell. Then we may move forward and solve a self-consistent dispersion 
relation. It is a well-known fact that in the high-temperature limit this is the case. 
Self consistent dispersion-relations have been solved in the high temperature limit 
without^, and also with a magnetic field using the Schwinger proper time method& 1 
The result, rich in physics, contains for example particle as well as hole solutions. 
With the general self-energy matrix presented here it is possible to perform the same 
analysis also in the Furry-picture. 

On the other hand the Furry-picture is particularly suitable in the strong field 
limit. In this limit the dominant contribution will come from intermediate electron 
states in the lowest Landau level (n = 0) only. In the limit {eB ^ ^ m^,p^,/i^} 
we find in the lowest Landau level, using only n = in the electron propagator Eq. (|^) 

Sr - a+mf In^ f ^ + A - l) , (13) 



TT 



4tt V 

a 



W( 1^2 In fill,/ 



(14) 



t2 - -a+m—\n(-] , (15) 



vr vm 



where we have defined m A = E — m — pj, as the deviation from the on-shell 
energy eigenvalue. The self-energy Sq = S™*^ + Sq^*^ + Sq^ has been split into its 



contributions from the vacuum, thermal fermions, and thermal photons, respectively. 
Similarly the gauge-fixing dependent parts of the off-shell self-energy are in this limit 



^vac(g) . . , f2eB\ 



oc fmaAln — —] , (16) 



Eo^'"(«) oc imaA\n[-] . {li 



T 



m 



Keeping only terms 0[ln^(2e5/m2)] and 0[\n{2eB/m'^) ln(T/m)], we may thus ne- 
glect the gauge fixing dependence. Also notice that all of the three different contribu- 
tions to the ^ dependent part of the self-energy are proportional to A. We have thus 
explicitly shown that the on-shell (i.e. A = 0) self-energy is gauge-fixing independent 
in this limit. The self-consistent dispersion relation obtained from Eq.(|ll]) reads in 
this limit 

(E70 - ^^7^ - m)a+Uo = dMncr+Uo , (19) 
where we have defined the "thermo-magnetic mass" 



SMb = ma/ (An) 



(20) 



The presence of 0"+ just signifies that we have already used that the wave-function in 
the lowest Landau level is proportional to the matrix Vo,py(x) = Iq^p^^x) diag[l, 0, 1, 0], 
according to Eq.(|). The above Eq.(^) is thus equivalent to the tree-level Dirac 
equation with the mass M = m + 5M.B- Unlike in the high temperature limit, this 
self-consistent dispersion relation thus only results in an energy shift 



E^^^ = J{m + 5MBy+pl ■ (21) 



To the lowest order in 5J^ b this is the same as would have been obtained from Eg . ([T^) . 
How large can then 5M.B be? In order for 5M.B ^ m we must have eB/m? k, 10^''', 
i.e. B K, 10^''T, an immensely large field. However, this makes our negligence of 
the ^ dependent parts more accurate then it appeared at first when keeping only the 
leading logarithms. The gauge fixing dependent part of the self-energy is proportional 
to A 6A4Bm/ Eq, that is small for most magnetic field strengths. 



5 Weaker Magnetic Fields 

In the case of an arbitrary magnetic field the sum over all Landau levels in the elec- 
tron propagator has to be considered. This problem has recently been solved^. Let 



us consider an on-shell electron in the lowest Landau level with vanishing momen- 
tum. We may now perform an expansion to linear order in e-B/A^, where A must be 
some energy scale, initially not known in this naive expansion. The result for the 
contributions to the self-energy from vacuum and thermal photons are, respectively, 

eB a 



AE2 



2m 2tt 
air 



eB 21X01 



m- 



3 VP? ^ 2m 9 

The thermal electron contribution is more involved. In the absence of a magnetic 
field we find the well-known result 



m^ 



(22) 
(23) 



[B = 0,p, = 0) = Am" 



Q/ roo I 

— / dkoQ{kl - m^)fF{ko)Jki - 



ko — 2m 

m{kQ — m) 



(24) 



To linear order in eB we find^ 



Am'^ 



+ 



eB a 
2m 3tt 



duo 



X 



2uj'^ + 2muj 



m 



- m^ 

m 



m^ 



+ 



eB a 
2m Sir 



duj 



{uo + m) 
2{jj^ — 3m'^u 



f+{uo) + 2m 



duj 



+ 



2m^ 



y/io"^ - m^ m'^{uj + m) 



-n 



(25) 



where some integrations by parts has been performed in order to make the result infra- 
red finite. The spin energy-shift for a particle of charge — e, with spin (s) magnetic 
moment /x = —e/{2m)g{s) in a magnetic field B = —Bz is AE = — /x ■ B. We may 
thus write 

AE = Am^'^-^^4 , (26) 



2m 2 ' 

where Am^'^ is the thermal mass, and 5g/2 is the anomaly of the magnetic mo- 
ment. Usually 6g/2 is obtained from the vertex correction (i.e. the triangle diagram). 
However, in the case of the contribution from thermal electrons the results differ. 
The discrepancy is caused by the external states being Landau levels, and not plane 
waves, as assumed when considering the vertex correction. Notice that the anomaly 
is obtained from the transverse part of the vertex correction, not included in the 
Ward-identity relating the vertex to the self-energy. 

Let us now consider the limit of vanishing temperature, but finite density of elec- 
trons. We may then perform the energy integral in Eq.(E3) and Eq.(|23) to find 



2m Sir m 

This may become very large at high densities, hke for example in the case of a neutron 
star, where /x/m ~ 10^, and T/m ~ 1. Comparing with the high temperature case, 
we see that the leading terms m{e^/mY has a similar origin to the corresponding 
m{eT /mY terms. Thus this contribution has been obtained from a "hard dense loop" , 
and higher loop corrections will only give sub-dominant contributionJil in powers of 
(e^/i/m). 



m 



2m^ 



m 



m 



(27) 



6 De Haas — van Alphen Oscillations 

Notice that the self-energy given in Eq.(P7D will become divergent as /i — >■ m"*" , T = 0. 
This is not physically acceptable, since the dense and thermal contribution should 
vanish in this limit. It seems as if our naive expansion in this case has been in powers 
of eB/lfi^ — m^), and thus is not valid for small densities. On the other hand, if 
2eB ^ fi^ ~ m? only a few lower Landau levels will be occupied by electrons. We 
thus have a situation very similar to the strong field limit, well suited for the explicit 
spectral decomposition. Split the thermal electron part of the self-energy into its 
contributions from the different Landau levels: 



N 
n=0 



e (n) 




Defining 



e (n) 



m[uj 



a 



m)/2eB — n, and G, 



n J y/m^+2eBn ^/up- 



X 



2eBn 

n-2 

(-a„)"-ie'^"Ei(a„ - ie) + 0^-2 



N = lnt[{fi^ -m^)/2eB] . {21 
1 , n > 0; 6r, = , n < —1, we find 



m ^ uj — m 



n\ {n 
'(n-2-/)! 



r 



.(29) 



In the limit /i — m only the lowest Landau level will contribute. Using the series 
expansion of the exponential integral, we find for /i^ — m? <^ 2eB 



a 



m— < In ( — - In 



TT 



eB\ f fJ' + V f^"^ — rn?' 



m- 



m 



H — m 



m 



2 - In 



m 



m 



(30) 

We can thus see that the electron contribution to the self-energy is vanishing as 
/i 171^ , T = 0, and furthermore that it is non-analytical in B. As (//^ — m'^) /2eB 
is becoming larger, consecutive Landau levels will cross the Fermi- surface, and start 
contributing. The contribution from the n-th Landau level is vanishing as /i — > 



\/m^ + 2eBn, so the self-energy is continuous. However, the derivative with respect 
to the chemical potential is diverging as — > y/m? + 2eBn. These sharp cusps will be 
smoothed out at finite temperature. There are thus oscillations in the self-energy at 
low temperatures, and Fermi-momentum squared of the order of the magnetic field. 
As a function of 1/eB, they have a quasi-periodicity of 2/{fi'^ — m^) = 2tt/A, where 
A is the area of an extremal cross section of the Fermi-surface. Such oscillations are 
well-known as de Haas-van Alphen oscillations in the magnetization of a Fermi-gas, 
but not earlier found in the fermion self-energy. 

For n > 1, a„ may become negative. This will cause an imaginary part in the 
self-energy according to Ei(— x — ie) = — Ei(x) + in, for x > 0. Let Te-^e--y denote 
the total decay rate for an electron in the Fermi sea decaying into an electron in the 
lowest Landau level, that not is supposed to be occupied. Then we find 

T,-^,-^ = 2lmAEf . (31) 

An ostensible contradiction appears here. The lowest Landau level, always below 
the Fermi energy, is supposed not to be occupied. However, perturbation theory 
rests on the assumption on adiabatic turn on/off of the interaction as t — > =fcxd. 
The external state is thus supposed to be separated from the heat and charge bath 
before and after the interaction is taking place, and not affected by the occupation 
of states in the medium. In Fig. |l] we show the real and imaginary parts of the dense 
electron contribution to the self-energy as a function of the inverse magnetic field. 
Unfortunately the cusps in the real part, as Landau levels with n > 2 are crossing 
the Fermi-surface cannot be distinguished on this scale. The cusp in the imaginary 
part at /i = m -|- 2eB/2m follows from Eg . (|29|) , but has no counterpart in ordinary 
de Haas-van Alphen oscillations. 
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Figure 1: The real part and the imaginary part (dotted hne) of the electron self-energy 
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